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Abstract 

The implications of shearfree condition on instability range of anisotropic 
fluid in f{R, T) are studied in this manuscript. A viable /(i?, T) model 
is chosen to arrive at stability criterion, where R is Ricci scalar and T 
is the trace of energy momentum tensor. The evolution of spherical 
star is explored by employing perturbation scheme on modified field 
equations and contracted Bianchi identities in f{R,T). The effect of 
imposed shearfree condition on collapse equation and adiabatic index 
T is studied in Newtonian and post-Newtonian regimes. 
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1 Introduction 


In a recent work [T], we have studied the effect of anisotropic fluid on dy¬ 
namical instability of spherically symmetric collapsing star in f{R, T) theory. 
Herein, we plan to explore the instability range of anisotropic spherically 
symmetric stars, considering shearfree condition. The role of shear tensor in 
evolution of gravitating objects and consequences of shearfree condition has 
been studied extensively. Collins and Wainwright [2] studied the impact of 
shear on general relativistic cosmological and stellar models. Herrera et ah 
lai worked out the homology and shearfree conditions for dissipative and 
radiative gravitational evolution. 

The features of gravitational evolution and its final outcome is of great 
importance in view of General Relativity (GR) as well as in modified theories 
of gravity. Shearfree collapse accounting heat flow is discussed in 0, it 
is established that shear plays critical role in gravitational evolution and 
may lead to the formation of naked singularities 0. It is mentioned in 0 
that occurrence of shearing effects near collapsing stars detains the apparent 
horizon leading to formation of naked singularity. However vanishing shear 
give rise to the formation of apparent horizon and so the evolving cloud ends 
in a black hole (BH). Thus, relevance of shear tensor in structure formation 
and its impressions on dynamical instability range of self gravitating body is 
well motivated direction of study. 

Stars shine by consuming their nuclear fuel, continuous fuel consumption 
cause imbalance between inward acting gravitational pull and outward drawn 
pressure giving rise to collapse [7] . The outcome of gravitational evolution is 
size dependent as well as other physical aspects 0 such as isotopy, anisotropy, 
shear, radiation, dissipation etc. In comparison to the stars of mass around 
one solar masses, massive stars tend to lose nuclear fuel more rapidly and so 
more unstable. Pressure to density ratio naming adiabatic index, denoted by 
T is utile in estimation of stability/instability range of stars. Ghandrasekhar 
0 explored the instability range of spherical stars in terms of T. 

Herrera and his contemporaries m-ng contributed maj orally in address¬ 
ing the instability problem in general relativity (GR), accompanying various 
situations, i.e., isotopy, anisotropy, shearfree condition, radiation, dissipa¬ 
tion, expansionfree condition, shearing expansionfree fluids. In order to 
achieve the more precise and generic description of universe, the dark en¬ 
ergy components are incorporated by introducing modified theories of grav¬ 
ity. Modified theories are significant in advancements towards accelerated 
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expansion of the nniverse and to present corrections to GR on large scales. 
The modihcations are introdnced in Einstein Hilbert (EH) action by indncing 
minimal or non-minimal conpling of matter and geometry im. 

Dynamical analysis of self gravitating sonrces in modihed theories of grav¬ 
ity has been discnssed extensively in recent years. The nnll dnst non-static 
exact solntions in f{R) gravity are stndied in [TB], Cembranos et ah [IB] 
stndied the evolntion of gravitating sonrces in the presence of dnst flnid. 
Instability range of spherically and axially symmetric anisotropic stars has 
been established in context of f{R) gravity [20112]) conclnding that devia¬ 
tions from spherical symmetry complicates the snbseqnent evolntion. 

Harko et ah [22] presented the f{R, T) theory of gravity as another al¬ 
ternate to GR and generalization of f{R) theory representing non-minimal 
matter to geometry conpling. The action in f{R,T) gravity inclndes arbi¬ 
trary fnnction of Ricci scalar R and trace of energy momentnm tensor T to 
take into acconnt the exotic matter. After the introdnction of f{R, T) grav¬ 
ity, its cosmological and thermodynamic implications were widely stndied 
[22]-[50] inclnding the energy conditions. Recently, we have stndied the evo¬ 
lntion of anisotropic gravitating sonrce with zero expansion [2T]. Herein, we 
are interested in exploration of shearfree condition implications on spherically 
symmetric gravitating sonrce in f{R,T) gravity. 

The modihed action in f{R, T) is as follows [22] 

+ ( 1 . 1 ) 

where C(m) denote matter Lagrangian, and g represents the metric tensor. 
Varions choices of £(m) can be taken into acconnt, each of which leads to 
a specihc form of hnid. Many people worked ont this problem in GR and 
modihed theories of gravity, stability of general relativistic dissipative axially 
symmetric and spherically symmetric with shearfree condition has been es¬ 
tablished in [321 [33] . Dynamical analysis of shearfree spherically symmetric 
sonrces in f{R) gravity is presented in [22] . 

The mannscript arrangement is: Section 2 comprises of modihed dynam¬ 
ical eqnations in f{R, T) gravity. Section 3 inclndes model nnder consider¬ 
ation, pertnrbation scheme and corresponding collapse eqnation alongwith 
shear-free condition in Newtonian and post Newtonian eras. Section 4 con¬ 
tains conclnding remarks followed by an appendix. 
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2 Dynamical Equations in /(i?, T) 

In order to study the implications of shearfree condition on evolution of 
spherically symmetric anisotropic sources, modihed held equations in f{R, T) 
gravity are formulated by varying action fll.ll) with the metric g^v Here, we 
have taken £(m) = p [27j, for this choice of £(m) modihed held equations in 
/(i?, T) gravity takes following form 

(/t + i)rir> - ps™/t + ^ 

+ (V^V^ — • (2-2) 

Here is energy momentum tensor for usual matter taken to be locally 
anisotropic. 

The three dimensional spherical boundary surface S is considered that 
constitutes two regions named as interior and exterior spacetimes. The line 
element for region inside the boundary E is 

ds^_ = r)dR — B^{t, r)dr'^ — r){d9'^ + sin^ 6d(j)^). (2.3) 

The line element for region beyond S is [1] 

ds\ = dv^ + 2drdu — r‘^{d6‘^ + sin^ Odcf)^), (2.4) 

where v is retarded time and M denote the total mass. 

The expression for anisotropic energy momentum tensor Tu^'^ is given by 

= {p + P±)VuVv - p±guv + {Pr - P±)XuXv, (2.5) 

where p is energy density, 14 describes four-velocity of the huid, Xu is radial 
four vector, pr and represent the radial and tangential pressure respec¬ 
tively. These physical quantities are linked as 

= = X^Xu = -1- (2.6) 

The shear tensor denoted by is dehned as 

^uv h(n;u) ®(tiI4) ^ 414 ), 
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(2.7) 






where is four acceleration and 0 is expansion scalar, given by 


an = V^u;v)V\ 0 = ( 2 . 8 ) 

Components of shear tensor are found by variation of Eq. fl2.7p and these 
are used to hnd expression for shear scalar in the following form 


1 





(2.9) 


where dot and prime indicate time and radial derivatives respectively. From 

shearfree condition we arrive at vanishing shear scalar, i.e., cr = 0, implying 
R = c 

B c 

It is worth mentioning here that the expansion scalar and a scalar func¬ 
tion described in terms of Weyl tensor and the anisotropy of pressure con¬ 
trols the departure from shearfree condition. Such function is related to the 
Tolman mass and appears in a natural way in the orthogonal splitting of 
the Reimann tensor [35]. It is obvious that pressure anisotropy and den¬ 
sity inhomogeneities have extensive implications on stability of the shearfree 
condition but it is not intuitively clear that their specihc combination affects 
stability [33|. Generically shearfree condition remains unstable against the 
presence of pressure anisotropy. Alternatively, one can consider such a case 
that pressure anisotropy and density inhomogeneity are present in a way that 
the scalar function appearing in orthogonal splitting of Reimann tensor van¬ 
ishes, implying non-homogeneous anisotropic stable shear free flow. Since we 
are dealing with fluid evolving under shearfree condition, so we shall make 
use of this condition while evaluating the components of held equations and 
also in conservation equations. 


5 


The components of modified Einstein tensor are 
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(2.13) 


The dynamical equations extracted from the conservation laws are vital in 
the study of stellar evolution. The conservation of full held equations is con¬ 
sidered to incorporate the non-vanishing divergence terms, Bianchi identities 
are 


GJK = 0, G^^xu = 0, (2.14) 

on simplihcation of Eq. fl2.14p . we have dynamical equations as follows 

p-P^ + [^ +fT]{^P + Pr + ‘2p±)^ +Zi{r,t) = 0, (2.15) 

r A' f A' C f' \ GA 

{P + Pr)fT + (1 + k) \p'r + 2^^ C I 

+/r(p'-£)+^ 2 (r,f)= 0 , (2.16) 

where Zi{r,t) and Z 2 {r,t) are provided in appendix as Eqs. fl5.1l) and (15.21) 
respectively. Deviations from equilibrium in conservation equations with 
the time transition leads to the stellar evolution, perturbation approach is 
devised to the estimate the instability range. 
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3 Perturbation Scheme and Shearfree Condi¬ 
tion 


We consider a particular f{R, T) model of the form 

f{R,T) = R + aR^ + \T, (3.17) 

where a and A can be any positive constants. Perturbation approach is 
utilized to estimate the instability range of spherical star with shear free 
condition. This scheme is utile in determination of more generic analytical 
constraints on collapse equation, rather to establish dynamical analysis of 
special cases numerically. Also, held equations are highly nonlinear differen¬ 
tial equations, in such scenario application of perturbation lead to benehcial 
observations. 

It is assumed that initially all quantities are independent of time and 
with the passage of time perturbed form depends on both time and radial 
coordinates. Taking 0 < e 1, the physical quantities and their perturbed 
form can be arranged as 


A{t,r) 

= Ao{r) + eD{t)a{r), 

(3.18) 

B{t,r) 

= Bo{r) + eD{t)b{r), 

(3.19) 

C{t,r) 

= Coir)+eD{t)c{r), 

(3.20) 

p{t,r) 

= Po{r) + ep{t,r), 

(3.21) 

Pr{t,r) 

= Pro{r) + epr{t,r), 

(3.22) 

p±{t,r) 

= P±oir) + ep~±{t,r), 

(3.23) 

m{t, r) 

= mo{r) + efh{t,r), 

(3.24) 

R{t,r) 

= Ro{r) + eDi{t)ei{r), 

(3.25) 

T{t,r) 

= To{r) + eD2{t)e2ir), 

(3.26) 

f{R,T) 

— [-Ro(^) + o:Roi'>') + -^^o] + ^iDi{t)ei{r)[l 



+ 2aRo{r)] + D 2 {t)e 2 {r)), 

(3.27) 

fn 

= 1 + 2aRo{r) + e2aDi{t)ei{r), 

(3.28) 

fr 

= A. 

(3.29) 


Considering Schwarzschild coordinate Cq = r and implementing perturbation 
scheme on vanishing shear scalar implies 


b c 
Bo r' 


(3.30) 
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Using Eqs. fl3.18p - fl3.29p and fl3.30p in dynamical equations i.e., Eqs. fl2.15p 
and fl2.16p leads to following expressions 


P + 


+ Ai-(2po + Pro + 4p_Lo) + y Zip 
Y r 


A' 


Al { Pr + +Pr [-f- ^ - 

Ao V^o r 


+2e^ - rfli I + B'^{Y) 


2 a Rq 
Y 


D = 0, 
2p-±\ 


(3.31) 
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(3.32) 


where Zip and Z 2 p are given in appendix, for the sake of simplicity we put 
Y in place of 1 -|- 2aRo and Ai = A -|- 1, assuming that Di = D 2 = D and 
Cl = 62 = e. Above mentioned perturbed dynamical equations and perturbed 
field equations shall be used to arrive at perturbed physical quantities such 
as p,pr and p±. 

The expression for p can be found from Eq. fl3.3ip . as follows. 


P 


^ 7-(3po + Pro + 4p_Lo) + Y Zip 

Y r 


D. 


(3.33) 


The Harrison-Wheeler type equation of state relate p and Pr, given by 


■p PrO - 

Pr = T -^- p. 

PO + PrO 

Putting p from Eq. fl3.33l) in Eq. fl3.34p . we found 

p^ = -V — ^+Xi-{3po + pro + 4:p±o) + YZip 
Po + Pro I Y r 


(3.34) 


D. 


(3.35) 


Perturbed form of field equation fl2.13l) yields expression for p±, that turns 
out to be 


Yc D Xp 

Px = {p-2a4^-- + 


A 

P±o ~ ^Po 


2ae Zo 
Y ^ Xi 


D. 


(3.36) 









































Zj, is effective part of the field equation given in appendix as Eq. fl5.5p . 

Substitution of andpd. from Eqs. fl3.33p . fl3.35p and fl3.36p in Eq. fl3.32p 
leads to collapse equation as under 
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f 2epo 




2p 


= 0 . 


(3.37) 


Matching conditions at boundary surface together with perturbed form 
of Eq. fl2.13p can be written in the simplihed form as follows 


D{t) - Z^{r)D{t) = 0, 

provided that 

Fc — 2Q;er 

The valid solution of Eq. fl3.38p turns out to be 

D{t) = -e^\ 


n , ^ _1_ \ 


(3.38) 


(3.39) 


(3.40) 


The terms of Z 4 must be constrained in a way that all terms maintain posi¬ 
tivity. Impact of shear-free condition on dynamical instability of N and pN 
regimes is covered in following subsections. 


3.1 Newtonian Regime 

In order to establish instability range in Newtonian era, we set po 3> PrOi 
Po 3> p±o and Aq = 1, Bq = 1. Insertion of these assumptions and Eq. fl3.40l) 
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in Eq. fl3.37p leads to the instability condition, relating usual matter and dark 
source contribution as under 


r < 



(3.41) 


where 



The quantities ^ip( 7 v) and Z 2 p(N) are terms of Zip and Z 2 p belonging to Newto¬ 
nian era. The gravitating source remains stable in Newtonian approximation 
until the inequality for T satishes, for which following constraints must be 
accomplished. 



The case when a —)■ 0 and A —?■ 0 leads to GR corrections and results for 
/(i?) can be retrieved by setting A —)■ 0. 

3.2 Post Newtonian Regime 

We assume Aq = 1 — ^ and Bq = 1 + ^ to evaluate stability condition 
in pN regime. On substitution of these assumptions in Eq. fl3.37p . we have 
following inequality for T to be fulhlled for stability range. 

Z4X5 + Xq + \pq{Xj + Y .Z'ip(pAr))q -|- XgXj — :;^Z^p]\f) + YZ2p{pN) 



(3.42) 
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Zip(PN) and Z 2 p[pN) are terms of Zip and Z 2 p that he in post-Newtonian era. 
The above ineqnality fl3.42p holds for positive dehnite terms and describe the 
stability range of snbseqnent evolution. The positivity of each term appearing 
in fl3.42p leads to following restrictions 
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4 Concluding Remarks 

In this manuscript, we carried out a study on implications of shear-free 
condition on stability of spherically symmetric anisotropic stars in f{R,T). 
Our exploration regarding viability of the f{R, T) model reveals that the 
selection of f{R,T) model for dynamical analysis is constrained to the form 
f{R, T) = f{R) + XT, where A is arbitrary positive constant. The restriction 
on f{R, T) form originates from the complexities of non-linear terms of trace 
in analytical formulation of field equations. The model under consideration is 
of the form f{R, T) = R + aFY + XT, representing a viable substitute to dark 
source and the exotic matter satisfying both viability criterion (positivity of 
radial derivatives upto second order). 

In f{R, T) gravity, the non-minimal matter geometry coupling include 
the terms of trace T in action fll.ip that is beneficent in the description 
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of quantum effects or so-called exotic matter. The components of modified 
field equations together with the implementation of shear-free condition are 
developed in section 2. Further conservation laws are considered to arrive 
at dynamical equations by means of Bianchi identities. These equations are 
utilized to estimate the variations in gravitating system with the passage of 
time. 

The complexities of more generic analytical held equations are dealt by 
using linear perturbation of physical quantities. Perturbation scheme induce 
signihcant ease in the description of dynamical system, rather to present sta¬ 
bility analysis by means of numerical simulations. The analytic approach we 
have employed here is more general and substantially important in explo¬ 
rations regarding structure formation. The perturbed shear-free condition 
together with the dynamical and held equations lead to the evolution equa¬ 
tion, relating T with the usual and dark source terms. It is found that the 
induction of trace of energy momentum tenor in action fll.ip contributes 
positive addition to T, that slows down subsequent evolution considerably. 

The outcome of gravitational evolution is size dependent as well as other 
physical aspects such as isotopy, anisotropy, shear, radiation, dissipation etc. 
The instability range for N and pN approximations is considered that im¬ 
pose some restrictions on the physical variables. It is observed that the 
terms appearing in T are less constrained for both the regimes (N and pN) in 
comparison to the anisotropic sources pQ. Thus, shear-free condition benehts 
in more stable anisotropic conhgurations. Corrections to GR and f{R) estab¬ 
lishments can be made by setting a ^ 0, A —)■ 0 and A —)■ 0 respectively. The 
local isotropy of the model can be settled by assuming Pr = p± = P- The 
extension of this work for shearing expansion free evolution of anisotropic 
spherical and cylindrical sources is in process. 
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5 Appendix 
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